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ABSTRACT 



We report on a 250 ks long X-ray observation of the supergiant fast X-ray transient 
(SFXT) IGR 116479^-514 performed with Suzaku in 2012 February. During this observa- 
tion, about 80% of the short orbital period (P or b'^3.32 days) was covered as continuously as 
possible for the first time. The source light curve displays variability of more than two orders 
of magnitude, starting with a very low emission state (10~ 13 erg cm~ 2 s _1 ; 1-10 keV) lasting 
the first 46 ks, consistent with being due to the X-ray eclipse by the supergiant companion. 

■ The transition to the uneclipsed X-ray emission is energy dependent. Outside the eclipse, the 

source spends most of the time at a level of 6-7 xl0~ 12 erg cm~ 2 s _1 punctuated by two 
structured faint flares with a duration of about 10 and 15 ks, respectively, reaching a peak 
flux of 3-4xl0~ n erg cm~ 2 s _1 , separated by about 0.2 in orbital phase. Remarkably, the 
first faint flare occurs at a similar orbital phase of the bright flares previously observed in the 
system. This indicates the presence of a phase-locked large scale structure in the supergiant 
wind, driving a higher accretion rate onto the compact object. The average X-ray spectrum is 
hard and highly absorbed, with a column density, Njj, of 10 23 cm~ 2 , clearly in excess of the 

■ interstellar absorption. There is no evidence for variability of the absorbing column density, 
\ except that during the eclipse, where a less absorbed X-ray spectrum is observed. A narrow 

Fe K Q emission line at 6.4 keV is viewed along the whole orbit, with an intensity which cor- 
relates with the continuum emission above 7 keV. The scattered component visible during 
the X-ray eclipse allowed us to directly probe the wind density at the orbital separation, re- 
sulting in p w =7x 10 -14 g cm~ 3 . Assuming a spherical geometry for the supergiant wind, the 
derived wind density translates into a ratio M^/oqo = 7xl0~ 17 M Q /km which, assuming 
terminal velocities in a large range 500-3000 km s -1 , implies an accretion luminosity two 
orders of magnitude higher than that observed. As a consequence, a mechanism should be at 
work reducing the mass accretion rate. Different possibilities are discussed. 
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1 INTRODUCTION 

IGR J16479^-514 is a hard X-ray transient discover ed by 
INTEGRAL on 2003, August 8-10 dMolkov et ai]|2003l) in the 
energy range 18-50 keV. Several X-ray flares were caught by 
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INTEGRAL/IBIS, displaying variable durations (from 0.5 to 
50 hr) and peak fluxe s (from 20 t o about 600 mCrab, 20-60 
keV;ISguera etalJl2005l , l2006t |2008l IWalter & Zurita Heras|[2007l 
iDucci et afl bOldl Recurrent outbursts were also observed by the 
Swift satellite fennea et alj|2005l [ Markward t & Krimmll2006l 
iRomano et ai]|2008llLa Parola et alj2009l.lBozzo et alj|2009h . 

The proposed counterpart (2MASS J16480656-45 12068, 
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IWalter et alj|2006l) wa s confirmed by an accurate localization ob- 
tained with Chandra dRatti et aT]|2010h and classified as a late O- 
type superg iant with a spectral type Q8.5I loca ted at a distance of 
4.9 kpc dRahoui et al.ll2008llChatv et alj|2008t) or a 09.5 lab star 
located at 2.81" kpc dNespoli et alJl2008l) ~The optical identifi- 
cation confirmed the initial classification of IGR J16479-4514 as 
a member of the new sub-class of high mass X-ray binaries, the 
supergiant fast X-ray transients (SFXTs), at first su ggested only 
based on the short duration of its hard X-ray activity dSguera et alj 
l2006h . 

IGRJ16479^514 is an eclipsing SFXT dBozzo et alj|2008bh 
and the on e with the narrow est orbit, showing an orbital period of 
3.32 days dJain et al.|[2009h . later refined to 3.3 193±0.0005 days 
dRomano et alj|2009h . 

X-ray long-term monitoring outside outbursts revealed that 
the source spends most of its time at a reduced level of X-ray 
emission (10 33 -10 34 erg s" 1 ), from 2 to 3 orders of magnitude less 
than the flare peaks, while for the remaining ~ 19% of the time, it 
is undetected, com patibly with being in eclipse ( Sidol i etalj|2008t 
iRomano et ai1l2009l) . 

The nature of the compact object is unclear (as in about a half 
of the members of the SFXTs class; for a recent review see Sidoli 
120101) but there is indirect evidence based on the spectral similarity 
with accreting pulsars suggesting that the X-ray source is a neutron 
star (NS). 

Here we report on the first X-ray observation which continu- 
ously (except that during the interruptions because of the satellite 
orbit) covers most of an orbital cycle of IGR J16479^45 14. Our 
main goal is to investigate the variability of the X-ray properties 
along a single orbital cycle. 



2 OBSERVATIONS AND DATA REDUCTION 



IGRJ16479-4514 was observed by Suzaku dMitsuda et al]|2007h 
between 2012 February 23 and 26. The observation (Obs ID 
406078010) was performed with the source loc ated at the X-ray 
Imaging Spectrometer (XIS; lKovama et alj2007l) nominal position. 
The XIS consists of four telescopes with a spatial resolution of 
about 2 arcmin coupled to four CCD cameras operating in the 0.2- 
12 keV energy range and with 18' x 18' field of view (1024 x 1024 
pixels; 1.05 arcsec pixel -1 ; Koyama et al. 2007). At the time of 
our observation only the two front illuminated (FI) XIS0 and XIS3 
(effective area: 330 cm 2 at 1.5 keV), and the back-illuminated (BI) 
XIS 1 (effective area: 370 cm 2 at 1 .5 keV) were operating. The three 
detectors were operated in Normal Mode with no window or burst 
option (all the pixels on the CCD are read out every 8 seconds). The 
other instrument aboa rd Suzaku is the Hard X-ray Detector (HXD; 
Takah ashi et al.l20 07). which consists of PIN silicon diodes (HXD- 
PIN; 10-70 keV) and Gd 2 Si0 5 Ce (GSO; 50-600 keV) scintilla- 
tors. The HXD was operated in the standard mode, with a time res- 
olution of 61 /is. The HXD is a collimated instrument and the PIN 
in particular has a 67.6' x 67.6' field of view (FWHM: 34' x 34'). 

Data reduction and analysis employed version 19 of the 
Suzaku Software included in the HEASOFT software package and 
followed the procedures described in the Suzaku ABC GuideQ 
Following standard practices, we excluded times within 436 s of 
Suzaku passing through the South Atlantic Anomaly and we also 
excluded the data when the line of sight was elevated above the 
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Earth's limb by less than 5°, or was less than 20° from the bright- 
Earth terminator. Moreover, we excluded time windows during 
which the spacecraft was passing through the low cut-off rigidity 
of below 6 GV. 

For the XIS we considered only events with GRADE = 0,2- 
4,6 and removed hot and flickering pixels using SISCLEAN; for the 
spectral analysis, the response matrices were generated with XIS- 
RMFGEN and using ray-tracing simulations with XISSIMARFGEN. 
With all the aforementioned data selection criteria applied, the re- 
sulting total effective exposure is roughly 136 ks for each XIS. The 
XIS events of IGR J 16479-45 14 were accumulated in each of the 
three XIS cameras within a circular region (3 arcmin radius) cen- 
tred on the target, while the backgrounds were estimated from an 
annulus with radii 5 and 7 arcmin. During the observation, the stan- 
dard criterion for the attitude determination requiring that the resid- 
uals are less than 0.005° could not be fullfilled, with residuals still 
remaining at about 0.008° for about 20 hours, in the temporal win- 
dow from February 25 01:00 through 22:00, but still smaller than 
the Suzaku point spread function (PSF) for usual point-like source 
analysis. 

In the HXD-PIN (no significant emission was detected in 
the GSO) the net exposure time, after dead-time correction (live 
time: 92.8%) is 141.9 ks. To subtract non-X-ray background (NXB) 
events that are included in t he HXD-PIN we use d the 'tuned' NXB 
synthesized background dFukazawa etai] [2009 ) . After subtracting 
the synthesized NXB from the HXD-PIN data, we also subtracted 
contributions of t he cosmic X- r ay ba ckground (CXB) using the 
model reported in iGruberetalJ Jl999h and of the Galactic Ridge 
X-ray emission (GRXE; IValinia & M arshall 199^). 

To ensure applicability of the x 2 statistics, the net XIS spec- 
tra were rebinned such that at least 20 counts per bin were present. 
All spectral uncertainties and upper-limits are given at 90% con- 
fidence for one interesting parameter. Data were analysed using 
FTOOLS version 6.11.1 and XSPEC version 12. In the spectral fitting 
we used the photoelectric ab sorption model PHA BS in XSPEC with 
the interstel lar abundances o f Wilms et al. ( 2000) and cross section 
table set of iBalucinska-Church & McCammon ( 1992). The three 
joint XIS0, XIS 1 and XIS3 spectra were fitted together, including 
constant factors to allow for normalization uncertainties between 
the instruments. For the timing analysis photon arrival times were 
corrected to the Solar System barycenter. 



3 ANALYSIS AND RESULTS 
3.1 Light curves 

The XIS barycentred and background-subtracted light curve of 
IGR 116479^-514 in the energy range 0.2-10 keV, is shown in 
Fig.[T] while the hardness ratio (HR) between net counts extracted 
in the hard (above 3 keV) and the soft (below 3 keV) energy bands 
is reported in Fig. [2] both versus time {upper panel) and in de- 
pendence of the source intensity (lower panel). Numbers in Fig.Q] 
indicate time intervals which show different source intensity be- 
haviour. At the beginning of the observation a very low intensity 
state is present (n. 1), consistent with being the X-ray eclipse pro- 
duced by the companion star. Indeed, foldin g the light curve on 
the refined orbital period P ort =286792±43 s dRomano et alj2009h 
and assuming the epoch 54547.05418 (MJD) as orbital phase (f>=0 
dBozzo etaL 1 20091) . we obtain the orbital phases reported as top 
x-axis in Fig.[T] The uncertainty on the epoch of 4>=0 can be esti- 
mated in about ^=±0.065, obtained extrapolating the error on the 
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Figure 1. Background-subtracted light curve of IGR J 16479^15 14 in the 0.2-10 keV energy range, obtained combining net source count rates observed by 
the three XIS units. In the upper x-axis, the orbital phase is indicated, assuming the orbital period P or j,=286792 s (Romano et al. 2009), assuming epoch 
54547.05418 MJD as orbital phase 0=0 (Bozzo et al. 2009). The bin size is 2048 s. Vertical dashed lines and numbers indicate the intervals of the eight 
time-selected spectra reported in Table|2] 



orbital period to the time of the Suzaku observation (432 orbital 
cycles between the epoch 54547.05418 MJD and the Suzaku ob- 
servations). 

The Suzaku observation did not cover the ingress time, so we 
are unable to determine the exact eclipse duration. We can at least 
constrain it to be between 46 ks (the duration of the eclipse at hard 
X-rays) and 143 ks (that is, 0.5 in phase, assuming a perfectly edge- 
on system with a tight orbit). An eclipse duration of about 0.6 days 
(52 ks) as proposed bv |jainetalj|2009l could imply that we might 
be seeing almost the full duration of the eclipse. After the eclipse 
(time interval n. 2) an intermediate level of emission is observed, 
lasting about 38 ks, during which the soft X-ray intensity below 3 
keV is indistinguishable from the emission during eclipse (Fig. [2)1, 
while at harder energies the source appears uneclipsed, although 
with a reduced intensity with respect to the uneclipsed emission 
displayed in the intervals n. 5 and n. 8. 

Two flares are present peaking at 2 counts s^ 1 and 
4 counts s _1 in the time intervals n. 4 and n. 7, respectively. In- 
tervals n. 3 and n. 6 can be considered the rise to these two flare 
peaks, while intervals n. 5 and n. 8 are characterized by a very vari- 
able, intra-flare emission, which covers about 40% of the orbital 
cycle. 



The X-ray light curve shows two sharp transitions from the 
eclipse towards the uneclipsed X-ray emission. To better estimate 
the time of the eclipse, we simply modeled the light curve with 
two ramp-and-step functions, characterized by seven parameters: 
the count rates of the three "steps" in the light curve and the four 
times, U, displayed in Fig. [3] The fit resulted in the following times 
measured from the beginning of the observation: ti=46 ± 1 ks, 
t 2 =47 ± 1 ks, t 3 =85 ± 2 ks, t 4 =88 ± 2 ks. The time t x corresponds 
to the date 55981.5131 ±0.0116 (MJD). 



3.2 Timing Analysis 

Barycentred light curves with 8 seconds bin time were produced for 
each unit of the X-ray Imaging Spectrometer (XIS0, XIS1, XIS3) 
in three energy bands (0.2-3, 3-10, 0.2-10 keV). We considered 
two time ranges: i) the entire duration of the Suzaku observa- 
tion (interval from 1 to 8 in Fig. QJ and ii) the part of the obser- 
vation excluding the eclipse and the transition to the uneclipsed 
emission (time intervals from 3 to 8 in Fig. [TJ. In order to search 
for spin periodicities in the XIS light curves, we used the Lomb- 
Scargle periodogram method by means of the fast implementation 
of Press & Rybicki (1989) and Scargle (1982), which is generally 
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Figure 2. Hardness ratio (hard/soft) variations with time (Upper panel) and with the source intensity {Lower panel). The bin size is 2048 s. Soft and hard 
energy ranges are below and above 3 keV, respectively. 



preferred for data set with gaps and unequal sampling. Periodici- 
ties were searched in the frequency range from i) 0.000013 Hz or 
ii) 0.000022 Hz (after which the sensitivity is reduced due to the 
finite length of the light curves) to 0.0625 Hz (corresponding to the 
Nyquist frequency of the data set). No significant and unambiguous 
evidence for coherent modulation was found in the periodograms. 



3.3 Spectroscopy 

We first analysed the joint XIS spectra extracted from the entire 
observation, corresponding to a net integration time of 136 ks. The 
net source count rates (0.2-10 keV) in the three XIS were the fol- 
lowing: 0.157±0.001 count s" 1 (XIS0), 0.144±0.001 count s" 1 
(XIS1) and 0.165±0.001 count s _1 (XIS3). Adopting an absorbed 
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Figure 3. Double ramp-and-step function fitting the initial part of the XIS 
light curve. The four times discussed in the Sect. 13. 1 l are indicated. 



power-law model resulted in a flat spectrum (photon index, T, 
of 1.35), a strong absorption (Nh~10 23 cm~ 2 ) and in positive 
residuals around 6.4 keV (Fig. [4j . Adding a narrow Gaussian line 
to the power-law continuum, we obtain a better fit, resulting in 
a line energy of 6.37 ± 0.03 keV, compatible with being pro- 
duced by Fe K Q fluorescence. The spectral parameters of the av- 
erage X-ray spectrum are listed in Table [TJ The average flux, F, 
of 1.16x 10 -11 erg cm -2 s" 1 (corrected for the absorption), trans- 
lates into an X-ray luminosity of 10 34 erg s _1 (assuming a distance 
of 2.8 kpc, iNespoli et al]|2008l) . Note howev er that, given the larg e 
uncertainty in the distance determination JNespoli et al.l [2008), 
this luminosity can range from 1.7xl0 33 erg s _1 (at 1.1 kpc) to 
8.2xl0 34 ergs" 1 (at 7.7 kpc). 

The hard X-ray spectrum extracted from the HXD-PIN field 
of view (FOV) is much brighter than the extrapolation at higher 
energies of the XIS power-law best-fit. Indeed, the HXD-PIN flux 
(~1.5xl0 -10 erg cm -2 s -1 , after taking into account CXB and 
GRXE) in the energy range 15-25 keV is ~15 times greater than 
that calculated extrapolating at high energies the power-law best 
fit to the XIS spectrum. Note however that this latter extrapo- 
lated flux should be considered as a conservative upper limit, since 
IGR J16479-4514 typically shows a spectral cut-off above 10 keV. 
This implies a strong contamination by other X-ray sources in 
the FOV, probably a residual contamination by the Z-track low- 
mass X-ray binary GX340+0, located at about 34 arcmin from 
IGR J16479-4514, and/or an unknown bright hard X-ray transient. 
For this reason, we will not discuss it further. 

The hardness ratio is variable (Fig. [2]l especially during 
eclipse, where the emission appears softer and/or less absorbed, so 
we next explored the variability of the spectral parameters in depen- 
dence of the orbital phase performing a time-selected spectroscopy 
extracting counts from the same intervals indicated in Fig.[T] Dur- 
ing the eclipse a remaining detectable flux is present, as typical in 
eclipsing high mass X-ray binaries, mainly produced by Thom- 
son scattering into the line of sight of the direct com ponent by 
free electrons residing in the supergiant wind (Haberl 1991). We fit 
the spectra with an absorbed power-law model, resulting in a less 
absorbed and softer spectrum during the eclipse compared to the 
other time intervals (Table [2} . The absorption during the eclipse 
is consistent with the interstellar reddening observed to the opti- 




Energy (keV) 

Figure 4. XIS spectrum extracted from the whole Suzaku observation. 
Count spectrum is shown (Upper panel), together with the residuals (in 
units of standard deviations) of the data to the absorbed power-law model 
(see Table[T] first column, for the spectral parameters). Lower panel shows 
the photon spectrum, graphically rebinned to better show the residuals 
around 6.4 keV. 



Table 1. Spectral results of the time averaged spectrum from the 
Suzaku/XIS observation. T is the power-law photon index. A Gaussian 
line is required at ~6.4 keV (see Fig. 0. Flux is in the 1-10 keV energy 
range in units of 10~ xl erg cm -2 s -1 and is corrected for the absorption, 
N H (in units of 10 22 
10~ 6 ph cm 



cm 



10 34 erg s" 



). The flux of the iron line, F; ine , is in units of 
2 s — 1 . Lx is the X-ray luminosity (1-10 keV), in units of 
for an assumed distance of 2.8 kpc. 



Parameter 


Power-law 


Power-law + Gaussian line 


Nn 


9.5 ±0.3 


9.5 ±0.3 


r 


1.33 ±0.04 


1.35 ± 0.04 


F, 

1 tine 




7.4 ±1.5 


EW (eV) 




70±10 


E Hne (keV) 




6.37 ±0.03 


Unabs. Flux 


1.16 


1.16 


Lx 


1.09 


1.09 


xS/dof 


1.018/2172 


0.991/2169 
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Figure 6. Fe K a emission line flux versus the flux in the energy range 7-10 
keV, corrected for the absorption. 



cal counterpart (Nespoli et al. 2008), while that measured outside 
the eclipse is significantly higher, and consistent, within the uncer- 
tainties, with being constant, at about Nh~10 23 cm~ 2 . The X-ray 
emission is harder when the source is brighter (X-ray flares), a typ- 
ical behaviour observed in SFXTs and in accreting pulsars. 

We next added a narrow Gaussian line to the single power-law 
model, resulting in the parameters reported in Tableland shown 
in Fig. [5] The resulting continuum power-law parameters are fully 
compatible, within the uncertainties, with those derived without the 
Gaussian line, so we do not report them again in Table [2] The en- 
ergy of the narrow line indicates K Q emission from neutral (or al- 
most neutral) iron. Its equivalent width (EW) has been calculated 
with respect to the power-law continuum, and its large value during 
the eclipse is due to the fact that it has been calculated with respect 
to the scattered component. The direct eclipsed X-ray emission is 
unknown, given the large source variability. The intensity of the flu- 
orescence Fe K a line is correlated with the unabsorbed hard X-ray 
flux (above 7 keV), at least outside the eclipse (Fig.[6]l. 

Finally, we adopted a power-law model together with a Gaus- 
sian line, but now fixing the power-law photon index to 1.35, the 
value obtained from the average spectrum. The fits are acceptable 
and the resulting values of the absorbing column density show now 
some variability, again with the less absorbed spectrum seen dur- 
ing the eclipse. More complex models are not required by the count 
statistics of the data and result in unconstrained spectral parame- 
ters. 



4 DISCUSSION 

The long Suzaku observation we have reported on here, allows 
us to investigate the properties of a SFXT as continuously as pos- 
sible within one single orbital cycle for the first time. This has 
been made possible by the short orbital period of IGR J 1 6479^45 1 4 
which is the SFXT with the shortest orbit known to date. The X-ray 
source shines with an average low luminosity Lx~10 34 erg s _1 , 
(if a distance of 2.8 kpc is assumed) and no bright flares (exceeding 
Lx~10 36 erg s _1 ) have been caught during this observation. Nev- 
ertheless, a high amplitude variability is observed, spanning more 
than two orders of magnitude (including the X-ray eclipse). The 



absence of signatures from highly ionizing illuminating continuum 
further confirms a low source X-ray luminosity. 

We can distinguish several types of variability on different 
timescales, especially evident in the energy range 3-10 keV. On 
long timescales of tens of ks we observe three main intensity states: 
the end of the eclipse by the supergiant companion (which we cov- 
ered for about 46 ks), a transition phase to the uneclipsed emission 
(time interval n.2 in Fig. 1, lasting 38 ks) and the intra-flare low lu- 
minosity emission which covers most of the remaining orbit (40%), 
lasting around 1 10 ks (time intervals n.5 and n.8). 

X-ray flares with different durations are present: two main 
low luminosity flares, lasting 10-15 ks, with a complex and struc- 
tured shape, punctuate the orbit, reaching peak luminosities of 4.4 
and 6.4xl0 34 erg s" 1 (at 2.8 kpc), separated by about 0.2 in or- 
bital phase. Their 0.2-10 keV spectra are harder than the intra-flare 
X-ray emission, as usually observed in SFXTs and accreting pul- 
sars. Even in the part of the orbit in-between the two main low 
luminosity flares, the source exhibits a large intensity variability 
on timescale of 1000 s. This low luminosity flaring activity had al- 
ready been unveiled in SFXTs (and in IGR J 16479-4 5 14) thanks 
to Sw ift/XKT long term monitoring observations (Sid oli et"al] 
l2008h . 

In the softer light curve the X-ray eclipse appears longer 
and there is no evidence for the presence of the transition phase 
(time interval n.2) clearly present in the hard light curve. Thus the 
eclipse egress appears to be energy depend ent. This behaviour is 
similar to that reported by I Jain et al. (2009) from the shape of the 
IGRJ16479-4514 folded light curves observed with RXTE and 
with SwiftfBAT, where the X-ray eclipse seems more evident at 
higher energies and with sharper transitions. There is also a sim- 
ilarity of the IGR J16479-4514 temporal variability near eclipse 
egress with what has be en observed in the HMXB 4U 1700-37 
(van der Meer et al. 2005 and references therein), showing a vari- 
able eclipse duration and a similar egress phase. The long duration 
(38 ks) of this phase roughly corresponds to a region as large as 
the companion radius (although a precise determination of both the 
donor radius and of the system inclination is lacking), probing the 
innermost structure of the supergiant wind. 

On shorter timescales of around 1000 s we observe two sharp 
transitions, indicated in Fig[3]by the four contacts U: (1) from the 
eclipse to time interval n. 2, which shows an intermediate level be- 
tween the eclipse remaining flux and the un-eclipsed flux, and (2) 
from time interval n. 2 to the first flare. Given that the compact ob- 
ject is point-like with respect to the supergiant companion, the first 
short transition likely probes the transition out off behind the stellar 
atmosphere, which, assuming an orbital velocity of a few hundreds 
km s" 1 , results in a linear size of around a few 10 10 cm. The second 
sharp transition is interestingly coincident with the rising phase to 
the first X-ray flare, suggesting a possible physical connection (see 
the discussion below). 

An iron line is detected along the orbit, with a line energy of 
6.4 keV, i ndicative of neutral iron or an ionization state lower than 
Fe XVIII dKallman & McCravll 19821) . A correlation is observed be- 
tween the iron fluorescent line flu x and the ill uminating unabsorbed 
flux above 7 keV, as expected foiouelll985h outside the eclipse 
(Fig. [6jl. The line flux observed during eclipse does not show ev- 
idence for a reduction, within the large uncertainties, with respect 
to the value detected during the intra-flare emission, so most of 
the fluorescing matter is located far away from the compact object. 
Note that we calculate the equivalent width of the line with respect 
to a single power-law continuum. If the fluorescing matter is opti- 
cally thin to the illuminating continuum, the luminosity of the iron 
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Table 2. Suzaku/XIS results of the time-selected spectroscopy (numbers mark the same time intervals displayed in Fig. 1) with an absorbed power-law 
model. Fluxes are in units of 10 -11 erg cm -2 s . Both observed (not corrected for the absorption) and unabsorbed (corrected for the absorption) fluxes are 
reported. Lx is the X-ray luminosity (1-10 keV), in units of 10 34 erg s —1 , for an assumed distance of 2.8 kpc. The absorbing column density, Njj, is in units 
of 10 22 cm" 2 . 
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Power-law 


Nn 

r 

Obs. Flux (1-10 keV) 
Unabs. Flux (1-10 keV) 
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The large EW during eclipse is calculated with respect to the scattered power-law continuum. 



K Q line depends on the spectrum and X-ray luminosity of the cen- 
tral source, on the c overing fraction o f the illuminated matter and 
iron abundance (e.g. lSako et alj fl999). Given the large uncertain- 
ties on the EW (especially huge during the eclipse phase) and the 
variability of the covering factor along the orbital phase, also con- 
sidering the likely presence of different large-scale structures in the 
supergiant wind which can be obvious, although unknown, sites of 
fluorescing material, it is not possible to gain further constraints on 
the properties of the reprocessing wind (also the supergiant surface 
becomes a reprocessor only visible at certain orbital phases). 

Considering the spectroscopy, a more physical description of 
the X-ray spectrum in an eclipsing HMXB involves the presence 
of two power-law continua, accounting for the direct and scattered 
component with the latter produced by Thomson scattering by free 
electron in th e companion wind and/or on the supergiant surface 
(Haber lll99ll) . Both the direct and the scattered components should 
have the same spectral slope, since Thomson scattering is energy 
independent. This implies that the only way to disentangle these 
two components is the different absorptions, the scattered compo- 
nent being less absorbed than the direct one. During the X-ray 
eclipse, where the direct component is unseen, the X-ray emis- 
sion should only be due to scattering. Indeed, fitting this spectrum 
with a power-law with a photon index fixed at the average slope 
of r=1.35, the eclipse shows the lowest absorption, consistent with 
the interstellar value ( iNespoli et alj2008h . Unfortunately, given the 
high absorption (local and interstellar) and the low statistics, we 
could not find better contraints to the spectral parameters adopting 
more complex models. 

Also the egress X-ray emission (n. 2) in IGR J16479-4514 
is likely dominated by the scattered component, given the high 



EW of the iron line compared to that measured during the out- 
of-eclipse emission (time intervals from n. 3 to n. 8). A possi- 
bility is that most of the direct X-ray emission in this orbital 
region is blocked by a large scale wind structure, as that in- 
voked to explain th e periodically recu rrent outbursts in the SFXT 
IGR Jl 1215-5952 dSidoli et al.ll2007l) . or by dense shells or other 
kind of inhomogeneities pr esent in hot massive sta rs winds (see, 
e.g.. lOskinova et alj|2012l or lLobel & Blommell2 008). This implies 
that most of the X-ray emission in this time interval is scattered 
X-ray emission, which is less absorbed because mainly produced 
by less dense wind material located off the orbital plane and far- 
ther away from the companion star. This could explain why the 
observed column density in time interval n. 2 is not larger than in 
the remaining (fully out-of-eclipse) part of the orbit. 

The likely presence of large scale structures in the supergiant 
wind is also suggested by the orbital phase of the first flare, which 
is compatible (within the uncertainties, see Sect. l3.lt with the lo- 
cation along the orbit of oth er bright flares obs erved in the past in 
IGRJ16479-4514 (<£~0.36: Ibozzq et alj[2009h . This suggests that 
the flares are triggered by a higher accretion rate during the pas- 
sage inside a phase-locked den ser wind component, like the one 
predicted bv lSidolietafl feofm This seems to be also indicated by 
the fact that, interestingly, the egress X-ray emission (n. 2) ends 
with the sharp transition to the first X-ray flare (see Sect. 13. j} , 
If this explanation is correct, after the first flare, the compact ob- 
ject lies in-between the observer and the dense wind structure, with 
the direct X-ray emission finally dominating the observed X-rays 
along the remaining part of the orbit. 

The second flare, spaced by about 0.2 in phase, could be rec- 
onciled by the crossing of a similar wind component, if the orbit is 
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Figure 5. Variability of the best-fit spectral parameters adopting an absorbed power-law model together with a Gaussian line at 6.4 keV (see Table |2). The 
observed (not corrected for the absorption) flux is in the energy range 1-10 keV and in units of erg cm - 2 s — 1 . The flux of the iron line is in units of 
photons cm -2 s — 1 . The orbital phase is indicated in the upper x-axis (see text). 



eccentric (as in the light curves simulated bv lDucci et ail2 009). or 
by the compact object approaching the periastron passage. Unfortu- 
nately, we lack any information about the other orbital parameters 
(e.g. eccentricity) to be able to confirm this hypothesis. 

The absorbing column density is always in excess of the inter- 
stellar reddening and does not show evidence for variability along 
the orbit (if also the slope of the power-law is mantained free), ex- 
cept that during the eclipse. On the other hand, if the power-law 



slope is fixed to the value of the average spectrum (r=1.35), vari- 
ability is present, with more absorption during the low luminosity 
flares. 

The only known orbital parameter is the period of 3.32 days. 
We did not observe the eclipse ingress, so we cannot constrain the 
companion radius, K op t- On the other hand, R opt can be estimated 
from the supergiant spectral type and requiring that the star does not 
overflow its Roche lobe. The spectral type of the supergiant (08.51 
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or 09.51) suggests a mass in the range M OP t =30-34 Mq and a 
radius of 22-23 R0 dMartins et alj|2003) . A further constraint on 
the stellar size comes from the f act that the comp anion should not 
overlflow its Roche lobe radius (Eggleton 1983). This suggests a 
donor mass around M op t=35 Mq and a stellar radius R op t=20 Rq . 
Adopting these values, the orbital period translates into an orbital 
separation of about 2.2 x 10 12 cm. 

We can use the X-ray eclipse to probe th e supergiant prop- 
erties, as follows (see, e.g., iLewis et alj|l992h . The average to- 
tal intensity during the eclipse is about 5% of the out-of-eclipse, 
intra-flare X-ray intensity, so the wind density can be obtained 
as n w = 0.05/(a<7T), where or is the Thomson scattering cross 
section and a is the orbital separation, which we assume as the 
characteristic path lenght through the system. This results in a 
wind density at the orbital separation, p w (a), of 7x 10 -14 g cm -3 . 
The mass continuity equation can be used to derive the ratio 
between the wind mass-loss rate and the terminal velocity, as 
M w /vao=A 7r a (a — R op t) pw(a), assuming a spherical ge- 
ometry for the outflow ing wind and a v elocity gradient for the 
wind velocity, f3, of 1 dCastor et alj|l975t) . This results in a ratio 
M u ,/-u 00 =7x 10 -17 MQ/km. For terminal velocities in the range 
500-3000 km s _1 , the mass loss rate is in the range M w =l- 
7x 10 MQ/yr. Assuming direct accretion from the wind, the ac- 
cretion rate can be estimated as M acc =(7rJ? 2 cc /47ro 2 )-M u) , where 
R acc is the accretion radius. This accretion rate translates into an 
X-ray luminosity Lx=3-15x 10 36 erg s _1 for the range of param- 
eters estimated before, which is two orders of magnitude greater 
than that we observe. 

Given the high wind density we have calculated, it is unlikely 
that the low luminosity in IGR J 16479^15 14 is due to the direct 
wind accretion. Our findings seems to agree with the recent re- 
sults by Oskinova et al. (2012), who have demonstrated that simple 
Bondi-Hoyle accretion from a clumpy wind over-predicts the ob- 
served X-ray variability in HMXBs. This suggests the presence of 
some mechanism able to damp the strong X-ray variati ons implied 
by the structured supergiant winds dOskinova et alj|2012l) . This can 
be due to the details (still poorly known) of the interaction of the 
accretion flow with the shocks in the accretion wake, leading, for 
example, to a transitional case of accretion regime, intermediate 
between Roche Lobe Overflow and direct accretion (as originally 
propos ed for the SFXT I GR J16418-4532, given its short orbital 
period dSidolietal.ll2012h . which displays a similar X-ray inten- 
sity variability). 

The role of the magnetospheric surface in reducing the accre- 
tion ra te onto the NS in SFXTs have been studied by several au - 
thors JBozzo et alJl2008j,|Pucci et alfeOlCtlshakura et alj|2012bl) . 
Bozzo et al. (2008) have invoked magnetar-like NS with slow pul- 
sations, to halt the accretion most of the time and allow only a 
residual flow of matter (pro ducing ^10 34 erg s - 1 ) by means of 
Kelvin-Helmholtz instabilitv. lShakura et al. (2012a) have proposed 
that the low X-ray luminosity state (~10 34 erg s" 1 ) sometimes 
observed in persistent accreting pulsars and in SFXTs could be due 
to subsonic quasi-spherical accretion onto slowly rotating pulsars, 
where the accretion of matter is mediated by a quasi-static shell 
above the NS magnetosphere. In this case, the accretion rate from 
the quasi-static shell depends on the ability of the plasma to enter 
the magnetosphere by means of plasma cooling. These authors sug- 
gest that transitions between low-luminosity states to bright flaring 
activity in SFXTs could be due to transitions between two different 
regimes of plasma cooling: from thermal plasma cooling to Comp- 
ton cooling dominated regime. Unfortunately, none of these mod- 
els can be confirmed to date, since both the pulsation period and 



the magnetic field in IGR J 1 6479^15 1 4 (and in most of SFXTs) are 
unknown. 



5 CONCLUSION 

The new Suzaku observations we have reported here allowed us 
to perform an in-depth investigation of the properties of the SFXT 
IGR J16479-4514 along one orbital cycle, and to obtain the follow- 
ing results: 

• the source spends most of the observation with an average 
out-of-eclipse (intra-flares) X-ray luminosity of 10 34 erg s _1 , with 
several kinds of variabilities on different timescales, but with no 
bright flares (exceeding 10 36 erg s _1 ) during the observed orbital 
cycle. Given long-term monitoring performed with previous mis- 
sions, this is very likely the most typical appearence of an orbital 
cycle in this SFXT; 

• the absorbing column density does not show evidence for 
variability, within the uncertainties, except from during the X-ray 
eclipse; 

• the remaining flux during the X-ray eclipse, produced by 
Thomson scattering, allowed us to estimate the wind density at the 
orbital separation, resulting in p w (o)=7x 10 -14 g cm -3 ; 

• assuming a circular orbit and a spherical geometry for the 
supergiant wind, the derived wind density translates into a ratio 
Mw/voo = 7x 10~ 17 MQ/km. Since the accretion luminosity im- 
plied by this ratio, assuming terminal velocities in the range 500- 
3000 km s~ , is at least two orders of magnitude higher than that 
observed, we can conclude that a mechanism mediating the accre- 
tion onto the putative neutron star in the system is likely to be at 
work to reduce the mass accretion rate. 
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